Mon. Not. R. Astron. Soc. OOO.mO('200x') Printed 2 February 2008 (MN I*TeX style file v2.2) 



CMB and foregrounds in WMAP first year data. 



G. Patanchon\ J.-F. Cardoso^, J. Delabrouille^, P. Vielva^ 

^ Department of Physics & Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, BC V6T IZl, Canada 
2 CNRS/ENST — 46, rue Barrault, 75634 Paris, France 

^ PCC — College de France, 11, place Marcelin Berthelot, F-75231 Paris, France 



Accepted ??. Received ??; in original form ?? 



ABSTRACT 

We perform a blind multi-component analysis of the WMAP 1 year foreground cleaned 
maps using SMICA (Spectral Matching Independent Component Analysis). We pro- 
vide a new estimate of the CMB power spectrum as well as the amplitude of the 
CMB anisotropics across frequency channels. We show that the CMB anisotropics are 
compatible with temperature fluctuations as expected from the standard paradigm. 
The analysis also allows us to identify and separate a weak residual galactic emission 
present significantly in the Q-band outside of the Kp2 mask limits, and mainly con- 
centrated at low galactic latitudes. We produce a map of this residual component by 
Wiener filtering using estimated parameters. The level of contamination of CMB data 
by this component is compatible with the WMAP team estimation of foreground resid- 
ual contamination. In addition, the multi-component analysis allows us to estimate 
jointly the power spectrum of unresolved point source emission. 

Key words: Cosmic microwave background - Cosmology: observations - Methods: 
data analysis 



1 INTRODUCTION 

\ The Cosmic Microwave Background (CMB) is one of the 
. most powerful probe of modern cosmology. The shape of 
' the spatial power spectrum of the small temperature fluc- 
tuations depend on the cosmological parameters describing, 
in the frame of the standard model, the mat ter content, the 
geom etry, and the evolution of the Universe jjungman et alJ 
Il996r . Since the first de tection of CMB an isotropies by the 
COBE sateUite in 1992 JSmoot et al.ll992^ . several ground- 
based and balloon-borne experiments have provided an ac- 
curate estimate of the power spectrum on a large range of 
angula r scales. The recent WMAP mission jBennett et alJ 
120033), after one year of data acquisition, provided measure- 
ments of the power spectrum with unprecedented accuracy. 

Measuring the CMB power spectrum is a difficult task. 
It requires a good characterization of noise contribution in 
the observations, as well as the subtraction of foreground 
astrophysical omission present at millimeter wavelengths 
iBouchet fc Gis pert 1999). Both noise and foreground emis- 
sions can significantly bias CMB power spectrum estimates 
if not well accounted for. 

The WMAP data is the best presently available obser- 
vation of the sky in the 20-90 GHz range. A careful separa- 
tion of CMB from foregrounds, however, is required in order 
to make the best out of these observations. In particular, 
isolating CMB from foregrounds is of extreme importance 



for measuring accurately the angular power spectrum C{£) 
of CMB anisotropies and for cosmological interpretation. 

The approach of the WMAP coUaboration consists 
in cleaning observations from foreground contaminations 
for performing a CMB power spectrum estimation on the 
cleaned maps. More specifically, galactic foregrounds are 
subtracted using templates obtained from external data sets. 
In addition, strong known extra-galactic point sources, as 
well as the region of the galactic plane, are masked prior 
to the analysis. The CMB spectrum is estimated using a 
weighted average of the cross-power spectra of cleaned ob- 
servation maps in the Q, V and W channels iHinshaw et alJ 
U0O33. This avoids biases due to detector noise, assuming 
that the noise of different detectors is uncorrelated. This es- 
timated spectrum is corrected for residual point source con- 
tamination by subtracting an estimate of the contribution of 
point sources in the cross power spectra. The level of resid- 
ual galactic contaminations is estimated by cross-correlating 
the maps with foreground templates. 

A multi-component approach of CMB 
spectral estimation has bee n proposed by 
IPclabroui lle. Cardoso. PatanchonI (|200^) . Additional 
publicati ons o n the subject give v ariants and details 
iCardoso et air 2002: P atanchon et alj [2003'l. The method, 
caUed SMICA (Spectral Matching Independent Component 
Analysis), is based on matching the cross- and auto-power 
spectra of the observed maps to a parametric model 
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described by the power spectra of all the astrophysical 
components, their relative amplitudes in the different 
channels, and the noise power spectra. It is a very flexible 
approach: depending on available knowledge, most of the 
parameters may be either estimated from the data or kept 
at fixed values. In particular, emission laws for all or some 
of the components can be estimated. 

There are several benefits to using SMICA in CMB 
analysis. 

First, SMICA is equivalent to a maximum likelihood es- 
timation of the CMB power spectrum if the observed maps 
are a linear mixture of Gaussian stationary components and 
noise. Therefore, if there are no foregrounds but only Gaus- 
sian CMB and noise, SMICA is expected to outperform 
quadratic estimates. 

Secondly, SMICA allows to estimate jointly the spa- 
tial power spectra and the amplitude of components in each 
channel (which are related to emission laws of components 
as well as calibration coefficients). In particular, the applica- 
tion on WMAP data allows us to check the following strong 
prediction of the standard model: CMB anisotropy should 
have a spectral emission law given, to first order, by a deriva- 
tive (w i th res pect to the temperature) of the blackbody law. 
iFixsenI l|200^ has measured the CMB anisotropy contribu- 
tion in the COBE FIRAS data by looking for the WMAP 
anisotropy template in FIRAS. The author shows the com- 
patibility of CMB anisotropy at large angular scales (> 5°) 
with temperature fluctuations in the Wien part of the spec- 
trum. Our blind estimation of the component emission laws 
provides us with a unique tool for extending this result to 
frequencies covered by WMAP and to all measured angular 
scales (including the first acoustic peak) with a remarkable 
precision as shown below. 

Finally, if there are foreground contributions in the ob- 
servations, SMICA is designed to detect them and allows 
their separation. It permits indeed to jointly estimate power 
spectra of multiple components in the data, to assess the 
significance of components, and eventually to separate the 
effect of all the emissions. 

In this paper, we investigate the existence of residual 
foreground emission in foreground-cleaned WMAP data. In 
particular, residual galactic emission resulting from an error 
in the subtraction of galactic templates can be present in the 
published maps. Such residuals may exist if external tem- 
plates of galactic emissions, extrapolated using a physical 
model to WMAP frequencies, are not quite representative 
of the actual galactic components. Multi-component analy- 
sis with SMICA allows to check for the existence of such a 
residual, and to quantify the level of contamination if any, 
without prior information on its power spectrum or on its 
amplitude across channels. 

In addition, the WMAP collaboration has shown that 
unresolved point sourc e emissions have a non-neKligible con- 
tribu tion to the data iHinshaw et al.lf2003l : iKomatsu et alJ 
[2003). This result has been i ndepen dently confirmed by 
iHuffenberger. Seliak. Makarovl (|2004 ). As will be shown, 
SMICA also yields a coarse estimate the level of point source 
emission in WMAP data. 

The paper is organized as follow: section |21 introduces 
the data used for the analysis, the simulations performed in 
order to check our results and the model of the WMAP ob- 



servations for SMICA. The method is described in section|21 
Section 01 presents the different results. Finally, conclusions 
are provided in section |3 



2 DATA 

2.1 WMAP data and input maps for SMICA 

The WMAP space probe, launched by NASA in 2001, is a 
large telescope for imaging the total emission of the sky at 
5 different wavelengths (or frequency channels) , with a reso- 
lution ranging from about 0.2 degrees to 0.9 d egrees (limited 
by diff raction), and with full sky coverage iBennett et alJ 

Eoosi). 

The data taken by WMAP has been made available to 
the scientific community after one year of proprietary pe- 
riod, and is freely available on a dedicated NASA CMB web 
site.^ The data consists in a set of 10 maps obtained by 
different detector pairs: four maps at 3.2 mm (W band at 
94 GHz), two at each 4.9 mm and 7.3 mm (V band at 61 
GHz and Q band at 41 GHz), and one at 9.1 mm and at 
13.0 mm (Ka band at 33 GHz and K band at 23 GHz). 
The data is pro vided in the HEALPix pixelli sation format 
of the sphere^ iCorski. Hivon. Wandelt]|l998l) . Foreground 
cleaned versions of the Q, V and W-band maps are also 
available. These maps have been used in the generation of 
the WMAP first-year CMB power spectrum by the WMAP 
team. The Galactic foreground signal, consisting of syn- 
chrotron, free- free, and dust emission, was removed using the 
3-band, 5-parameter template fitting method described in 
iBennett et all fe003bV Galactic templates come fro m exter- 
nal d ata sets: the 408 MHz synchrotron map (Haslam et alJ 
Il982h . t he predicted dust emission at 94 GH z using the FDS 
model frinkbeiner. Davis fc Schl egeJll99tt l. composite Ha 
map (F inkbcincr 2003) and the galactic reddening E(B-V) 
JSchlegel. Finkbeiner fc Davisiliggi) . 

For our analysis, we use the eight individual foreground 
cleaned maps. We partially correct the maps from beam 
smoothing effects, so that every map gets the effective reso- 
lution of a reference map (we choose W3). We simply multi- 
ply the coefficients of the spherical harmonic decomposition 
by the inverse of beam transfer function ratio. 

We consider two different sets of maps corresponding to 
different masking. For the first set (hereafter map set I), we 
apply the Kp2 galactic mask provided by the WMAP collab- 
oration, canceling about 15% of the pixels. For the second 
set (hereafter map set II), we mask the pixels corresponding 
to galactic latitudes lower than 40 degrees. For both sets, we 
mask the strongest point sources using the mask provided 
by the WMAP team. In all cases masks are apodised for a 
smooth transition between and 1 on a scale of 30 arcmin. 

The prior correction of beam smoothing effects permits 
the application of SMICA using directly the pseudo cross- 
and auto-power spectra of observation (computed from par- 
tially covered maps). The true component power spectra can 
be obtained from these pseudo power spectra afterwards. 



^ http://lambda.gsfc.nasa.gov/ 

^ http:/ /www. eso.org/science/healpix/ 
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2.2 Simulations 

In order to validate the method, assess error bars on the es- 
timated parameters, and check for systematic errors in the 
analysis, we have generated 100 full-sky simulations of the 
WMAP observed maps with 6.5 arcmin pixels (HEALPix 
nside = 512). Each set of simulations consists in 8 maps re- 
producing the observations in the Q, V and W bands. Maps 
contain synthetic CMB anisotropics degraded to the resolu- 
tion of the detectors, assuming a symmetric beam pattern 
and using the transfer function published by WMAP. They 
contain anisotropic white noise at the expected level. CMB 
anisotropies are generated using the CMBFAST software 
jZaldarriaea fc Seliakll2000h with current 'concordance' cos- 
mological parameters fep ergel et al ]l2003f). The simulated 
maps are partially deconvolved and the same mask as used 
for real data is applied. 

No measurable bias has been observed in simulations. 
Error bars on all parameters (CMB power spectrum as well 
as mixing parameters) obtained from the dispersion over 100 
independent simulations have been found to be in very good 
agreement (precision better than 10%) with the closed form 
expression of eq. (1131 . validating the use of the latter for 
final error bar estimates. 



2.3 Model of the WMAP observations 

The sky emission at WMAP frequencies is well described 
to first order by a linear superposition of the emissions of a 
few processes: CMB anisotropies, galactic foregrounds (syn- 
chrotron, dust, free-free. . . ), point source emissions and Sun- 
yaev Zel'dovich (SZ) effect. After subtraction of the galactic 
foreground templates and masking, the WMAP data be- 
tween 40 and 94 GHz can be modeled as noisy linear mix- 
tures of CMB anisotropies, unresolved point sources and 
possibly small residual galactic emissions (no significant SZ 
effect is expected to be pr esent in WMAP observations 
jBennett et alJ l2003b^ , and iHuffenberger. Seliak. Makarovl 
i200 4r) found no evidence of SZ effect using a multi-frequency 
approach). We assume that the emission laws of components 
are independent of the position [6, (j)) on the sky. Although 
this is only an approximation, it holds exactly for the CMB 
and for galactic residuals which are significant at one fre- 
quency only. Thus, observation i is modeled as 



Xi{e, 0) = ^ AicSc{e, 4)) + ni{e, 



(1) 



where Sc is the spatial distribution of component c, rii is the 
noise of observation i and Aic is the amplitude of component 
c in map i given by 



(2) 



where Wi is the spectral band of detector i and gc(i^) is 
the emission law of component c. In matrix- vector format, 
model Q reads 



x{e,(t)) = As{e,(f>) + n{e, 



(3) 



with an A'^,; x A*'c mixing matrix A. 

This simple model does not reflect the fact that reso- 
lution depends on the frequency band. In our analysis, this 



effect, as well as the impact of the masks, is accounted for 
in the spectral domain: see section 13.41 , for a sketch of the 
correction of beam and coverage effects. 

On the statistical side, we shall assume statistical inde- 
pendence between components, between the noise contami- 
nations of different detectors, and between the latter and all 
the components. 

The standard model of cosmology predicts that CMB 
anisotropies are small temperature fluctuations of a pure 
blackbody spectrum. A flrst order expansion around T = 
2.726 K gives: 



dT 



(4) 



where Bi,{T) is the Planck law for the emission of a black- 
body at temperature T. 

Whereas the assumption that Aic does not depend on 
the position is a good approximation to first order for galac- 
tic components, it is not the case for point source emissions. 
Nevertheless, since the brightest point sources are masked 
and since the frequency dependence of most of contributing 
radio sources belong to the flat population, their emission 
law can be roughly described as: 



(5) 



with /? « -2 llToffo latti et al.l ll998t) . Note that unresolved 
point sources in the WMAP data have alr eady been reported 
at a very low level Ikinshaw et alJl2003l) . 



2.4 Spectral statistics 

The spectra and cross-spectra at frequency I oi x{9, (j)) (con- 
sidered as an isotropic A'^^-dimensional random field on the 
sphere) is the Nd x Nd spectral matrix R{£) 



R{e) = {x{£,m)x{e,m)'' ) 



(6) 



where (■) denotes an ensemble average, superscript ^ denotes 
transposition and the x{l, m) axe the coefficients of the field 
on the basis of real spherical harmonics. This average value 
is independent of m because of isotropy. 

In the linear model of eq. J^J, independence between 
components and noise implies that 



R{t) = AC(t)A^ + N{1) 



(7) 



where C{£) and N(£) are the spectral matrices for the com- 
ponents and the noise respectively. They are both diagonal 
matrices as a consequence of the independence assumption 
between components and between noise contaminations. 

In practice, spectral matrices are estimated by aver- 
aging over frequency bins. Typically, one considers Q fre- 
quency bins with the q-ih bin (1 < q < Q) containing all 
harmonic modes {£,m) such that Iminiq) < t < ^max(f?). If 



{21 -\- 1) denotes the number of modes in 



the gth harmonic bin, then the empirical spectral matrix 

«max(<j) I 

Rq = - ^ ^ a;(£,m)x(£,m)t (8) 
is the natural estimate of the average spectral matrix 
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fmax(9) 



Tin 



^=*min(9) 

The latter inherits its structure from ^ as: 



Eg = ACgA' + iV, 



(9) 



(10) 



with average spectral matrices Cq and Nq related to C{£) 
and N{£) as in definition @. 



3 SPECTRAL ESTIMATION AND 

COMPONENT SEPARATION WITH SMICA 

The goal of our analysis is to identify and separate the com- 
ponents present in WMAP maps, to evaluate the amplitude 
of CMB anisotropics as a function of observation wavelength 
(as given by the corresponding column of A), and to provide 
accurate estimates for power spectra of the CMB and other 
components. This section briefly describes a multi-detector, 
multi-component analysis method. It is a maximum likeli- 
hood method based on a model of statistically independent 
components. Since it can be understood as a spectral match- 
ing technique, it is dubbed "SMICA", standing for "spectral 
matching independent co mponent analysis" . More details on 
SMIC A can be found in iDelabrouille. Cardoso. Patanchonl 
J2003h. 



3.1 Spectral matching 

The SMICA analysis technique consists in minimizing a 
measure of the mismatch between the empirical covariance 
matrices Rq and their model counterparts Rq. Minimization 
is conducted with respect to any relevant set of parameters 
describing the covariance matrices Rq . The maximal param- 
eter set is made of: the entries of the mixing matrix A, the 
average power [Cg]cc of the c-th component in the g-th bin, 
and the average power [Nq]ii of the noise on the i-th detector 
in the q-th frequency bin. 

One may choose to fix some of these quantities (prior 
knowledge) and to estimate the remaining ones, either freely 
or under some additional parametric constraints (see section 
14. H . In any case, denoting by S a set of free parameters defin- 
ing the values of Aic, [Cq]cc and [A^qjii, we obtain an estimate 
of 9 by minimizing a joint spectral mismatch, quantified by 
a weighted sum 



m 



9 = 1 



{r„r,{o)) 



(11) 



data do come from an isotropic random field. Also, the con- 
nection with the likelihood criterion suggests a simple opti- 
mization strategy using the EM (expectation-maximization) 
algorithm. In practice, minimization of the spectral mis- 
match mil l is achieved with the EM algorithm followed by 
few steps of a descent method (BFGS algorithm) to speed- 
up the final convergence. 



3.2 Error estimation 

The SMICA estimator being a maximum likelihood estima- 
tor, the asymptotic variance-covariance matrix of the esti- 
mates is given, when the model holds for some value 6o of 
the parameters, by 

Cov(e) = {{{6 - eo)ie - eo)*}> ~ J(eo)-' (i3) 

where the Fisher information matrix J{9q) is, in our model 



R^'^R^'^i (14) 



Any sensible measure D(-, ■) of matrix mismatch could be 
used to obtain a consistent estimate 6 of the 6 parameter as 
6 = argmiufl "I>(S). We adopt 

D{Ri,R2) = i [trace (jJiiia"') - log det(i?iiij') - iVd] (12) {nqD{Rq, Rq{9))) 

because then criterion 1)11^ is, up to a constant, equal to 
minus the log-likelihood of the data in a simple Gaus- 
sian is otropic model f see IDelabrouille. Cardoso. Patanchonl 
l)200.j) for a derivation). Hence, minimizing criterion l|ll|l 
is equivalent to maximizing the likelihood. This fact guar- 
antees good performance (at least for large nq) when the 



[J{9o)]ij = Tig trace 



with all quantities evaluated at point 9o. In practice, J{9o) is 
approximated by replacing, for each frequency bin q, matrix 
Rq by Rq and matrix ^^(9o) by -^{9). 

Departure from stationarity^ and Gaussianity does not 
introduce bias in parameter estimation but is likely to in- 
duce a larger estimation error than predicted by the above 
formula. The accuracy of the error bars prediction can be 
checked thanks to simulations. 

The galactic and point source mask has the effect, to 
first order, of decreasing the effective number Uq of modes 
in each bin by a factor equal to the fraction of sky coverage. 
This is accounted for by multiplying every parameter error 
estimate by the inverse square root of this factor. This pro- 
cedure has been validated with the help of numerous Monte- 
Carlo simulations. 



3.3 Mismatch measure 

The SMICA approach includes a built-in goodness of fit: 
the departure of the spectral statistics Rq from the best fit 
model of independent components is quantified by the di- 
vergence D{Rq, Rq{9)) which should be statistically small. 
For a known Nd x Nc mixing matrix, if the model of ob- 
servations is correct, then 2nqD{Rq, Rq{9)) is asymptoti- 
cally (i.e. for Uq large enough) distributed as a with 
Nd{Nd -f l)/2 - {Nc + Nd) degrees of freedom. This value is 
obtained by subtracting the number of adjustable parame- 
ters per frequency bin {Nc auto-spectra plus Nd noise levels) 
from the Nd{Nd + l)/2 degrees of freedom of a symmetric 
matrix. In particular: 



Nd{Nd + 1) 



{Nc + Nd) 



(15) 



In particular, as a consequence of WMAP scanning strategy, 
the noise variance per pixel on the observed maps is higher by a 
factor of the order of 4 around the ecliptic equator as compared 
to the pole regions. 
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3.4 Beam and coverage effects 

We discuss the corrections needed to account for beam and 
coverage effects. 

In first approximation, if a mask covers a fraction a of 
the sky, an effective number auq should be used in place of 
Uq. However, since masking also introduces mode coupling, 
a better approximation is desired. W e follow the MASTER 
formalism developed bv lHivon et alJ (l2002ll . The idea is as 
follows: let s{9, <j)) be an isotropic Gaussian random field 
with harmonic spectrum c{£), of which only a masked version 
s{6,<j)) ~ w[9,()))s{6,(f)) is observed. Denoting s{£,m) the 
harmonic coefficients of s{9, (j)), we have 

^ ^ (|s(Am)r> = ^M«,c(/) (16) 

m^-e V 

where the coefRcients Mui depend only on the mask (and 
thus can be computed independently of the data) . Hence, if 
M^i is known, so is the bias introduced by the mask on the 
harmonic spectra. 

In practice, we proceed as follows. In a first step, the 
input maps are masked, transformed to the spherical har- 
monic domain and brought to a common resolution (see sec- 
tion l2.1ll . Next, bin-averaged empirical spectral matrices are 
computed according to eq. JSJ and a model of independent 
components is adjusted to them by minimizing IIH . The 
resulting bin-averaged harmonic spectra of each component 
is then corrected by inverting the bin-averaged version of re- 
lation 1161 1. This last stage also incorporates the correction 
of the common beam pattern. 



3.5 Component map separation 

Ideally, component maps would be estimated by applying a 
Wiener filter W = [A^N'^A + C~^]~M*iV-i to the obser- 
vations. This solution maximizes the signal to noise ratio on 
each individual component map (here noise means detector 
noise plus other astrophysical component emissions). In the 
limiting case where noise is small as compared to compo- 
nent signals, is negligible and the Wiener filter yields 
unbiased (in the sense that WA = 7) estimates of the maps. 
In poor signal to noise regimes, the signal is attenuated to 
suppress noise contamination in the reconstructed maps. 

In practice, the Wiener filter is applied in the frequency 
domain, using estimated values A, N and C of the parame- 
ters. The harmonic coefficients of the estimated components 
are obtained as 

s{l,m) = [A'Nq^A + Cq'^r'^A'N-^x{l,m) (17) 
when £min(q) < t < ^max(a). 



4 APPLICATION & RESULTS 

We now present a SMICA analysis of WMAP data. The two 
data sets described in section r2.1l (map set I corresponding 
to sky regions outside of Kp2 mask, and map set II corre- 
sponding to galactic latitudes higher than 40 degrees, with 
the strongest point sources masked in both cases) are used 
for two independent analyses. 



4.1 Choice of parameters for SMICA 

In a preliminary analysis, we have looked for the total imm- 
ber of components required by the map set I without impos- 
ing constraints neither on the amplitude of the estimated 
components nor on their spatial power spectra. In this to- 
tally blind analysis, we found two significant components. 
The first component, dominant in all channels, is clearly 
identified as CMB. The second, weaker by several orders of 
magnitude and essentially significant in the Q band (it is de- 
tected to a lesser extent in the V and W bands) is thought 
to be a mixture of residual galactic emissions and unre- 
solved point sources. These last two processes cannot be re- 
liably separated without introducing additional constraints 
because of nearly proportional mixing matrix columns (both 
components dominate in the Q band and are almost negli- 
gible in the others). 

In order to differentiate between residual galactic emis- 
sion and point sources, we need to introduce some physical 
knowledge in the form of constraints on the system. There- 
fore, in most of the subsequent analysis, unless explicitly 
stated otherwise, we adjust a model with three components, 
one of them being strongly constrained to capture point 
source contributions as follows. The emission spectrum of 
point sources is well described, at WMAP frequencies, by 
a power law with a spectral index l3 ~ —2 (section 12.31 
and their spatial power spectrum is expected to be almost 
ffat since the effect of clustering is negligible in practice at 
radio frequencies |Gonzalcz-Nucvo ct al. 2004). Hence, for 
one component, meant to be residual point source emission, 
we fix the mixing parameters (the column of the mixing ma- 
trix) to Ai.ps ~ ivi/vo)^'^ in RJ temperature units, and we 
constrain the harmonic power spectrum to be fiat. Only one 
parameter, its amplitude, is left free to match the contribu- 
tion of this component. 

In summary, we match spectral matrices with three 
components. For two of them, meant to be CMB anisotropies 
and residual galactic emission, no constraint on parameters 
are enforced (both the mixing parameters and the power 
spectra are determined exclusively from the data). The last 
component is constrained to have the emission law and spec- 
tral shape of point sources; only its amplitude can be ad- 
justed (but see figure |7| and the related comments). Finally, 
regarding noise, the average noise power is freely estimated 
in each map and each frequency bin. 

4.2 CMB anisotropies 

Figure shows the estimated CMB power spectrum after 
correcting for partial coverage and beam and pixel transfer 
functions. It displays a peak around I — 200 (first acoustic 
peak) and a second peak around £ — 550, both compatible 
with the measurement announced by the WMAP team. 

SMICA and WMAP team power spectrum estimations 
show an excellent agreement for most multipoles. For multi- 
poles between £ = 2 and £ — 290 the difference between the 
two power spectra is much smaller than error bars (except 
for bins centered on ^ = 190 and £ = 210 where the difference 
is of the order of the error bars) , but is larger than statistical 
errors after removing cosmic variance contribution. This can 
be explained (at least partially) by the small differences in 
the sky coverage due to the apodisation procedure we apply 
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Figure 1. The CMB spatial power spectrum measured witli SMICA (in red) as compared to tlie publislied WMAP 1 year spectrum (in 
blue). The error bars from SMICA are computed from the Fisher information matrix using the parameters at convergence (and not the 
parameters of the model). This explains the small error bar on the first bin of the SMICA estimate. We also plot the difference of the 
power spectra for better comparison. The red curves correspond to ztlcr error bars on the power spectrum obtained with SMICA. 
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Figure 2. Measured amplitude of CMB fluctuations in each de- 
tector (in red) relative to temperature fluctuation units as given 
by the calibration on the dipole. The amplitude of the fluctuations 
are normalized to 1 on average. Data points at the same frequency 
have been slightly offset in abscissa for readability. Blue lines de- 
limit calibration 1 sigma errors provided by the WMAP team. 
The errors on estimated mixing parameters are statistical errors 
computed from the Fisher information matrix at convergence. In 
order to provide an error bar on every CMB mixing parameter, 
errors are not marginalized on CMB power spectrum estimate 
(there is a degeneracy between power spectrum amplitude and 
mixing parameters normalization). 
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Figure 4. Cross-correlation spectra between our CMB map with 
other available CMB maps from the WMAP data: the ILC map 
published by the WMAP team (Bennett et al. 2003b), the com- 
bined map used for non-Gaussianity studies (e.g., Komatsu et al. 
2003) and the two maps {cleaned and Wiener filter) obtained by 
Tegmark et al. (2003).) 



on the maps. For larger multipoles, the two estimates show 
somewhat larger differences which may be due, in part, to 
the two different weighting schemes used by the two meth- 
ods for e > 200. For SMICA, no weighting of the pixels 
depending on the noise variance per pixel is applied prior 
to the analysis. On the contrary, the WMAP team applies 
different weighting schemes for £ > 200 (a weighting pro- 
portional to 1/noise^ for £ > 450 and a transition weighting 
for 200 < £ < 450). Hence, discrepancies between the two 
estimates are expected to be larger for £ > 200 since the two 
data sets are not quite identical. This may account for the 
small observed differences. We note in the passing that the 
small dip at the top of the first acoustic peak in the WMAP 
team estimate has disappeared in the SMICA estimate. 

Figure |21 shows the estimated amplitude of CMB 
anisotropics for all the observation channels, as given by 
the CMB column of the estimated mixing matrix A. We 
find an emission law compatible with the expected deriva- 
tive of a blackbody, to excellent accuracy. A fit of the form 
gives a = -0.0067 ± 0.0058 compatible with at about 
1.15 sigma. Refined calibration and systematics testing, as 
well as the second year data, will reduce the error bar and 
clarify the significance of a possible deviation from the ex- 
pected blackbody derivative. Note that, given the error bars, 
the accuracy of the measurement of the CMB emission law is 
limited by detector calibration uncertainty rather than intrin- 
sic statistical errors! The present measurement shows that 
the emission law of the anisotropies is the same as that of 
the dipole, which is itself known (under the assumption that 
it is essentially due to a kinetic effect) to be the derivative 
with respect to temperature of the CMB blackbody emission 
measured with FIRAS. 



Figure El displays a reconstruction of the CMB 
anisotropy map by Wiener filtering (equation 1171 using the 
parameters estimated with SMICA. The now familiar CMB 
anisotropy pattern is clearly visible. 

We have compared our recovered CMB map with other 
available CMB maps from WMAP: the Internal Linear 
Combination (ILC) m ap published by the WMAP team 
jBennett et alJl2003bh . th e combined map used for non- 
Gaussianity analysis (e.g., iKomatsu et al.l 12003*)) and the 
two maps (cleaned an d Wiener filter) obtained by Tegmark 
jTeemark et alJl2003t) . 

As seen on figure |1] we find very high correlation 
(~ 99%) with the two maps proposed by Tegmark as well 
as with the map used for non-Gaussianity studies up to 
£ « 200. The correlation slightly decreases down to « 95% 
for multipoles larger than 200 (see Figure llj. The cross- 
correlation with the ILC map is also around of 99%, but 
only up to multipoles £ ~ 20. For smaller scales, the correla- 
tion with the ILC map decreases down to « 60% at ^ « 200 
and then it increases up to « 90% at the smallest scales. The 
fact that the correlation with the ILC is smaller than with 
the other CMB maps is not surprising since, as explained 
by the WMAP team, the ILC was developed for foreground 
analyses and it is far for being the best CMB reconstruction 
(Bennett et al. 2003b). 

We have also tested possible contamination caused by 
unsubtracted foreground emissions by cross-correlating our 
CMB reconstruction with templates of different Galactic 
emissions (synchrotron, thermal dust and free- free). We find 
no significant contamination, since the correlation level is be- 
low the dispersion of the casual correlations between those 
templates and CMB Gaussian simulations (5 — 10%). 
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Figure 3. Map of CMB fluctuations obtained after Wiener flltering of the data using parameters estimated with SMICA. 
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Table 1. Estimated relative amplitude of the second component 
in each of the observed map in Rayleigh-Jeans (RJ) temperature 
units and associated errors. The normalization of parameters is 
fixed with respect to Ql map. Error on parameters are correlated 
at the 20% level. 



4.3 Residual galactic emission 

TableQgives the emission law (the column of the mixing ma- 
trix) estimated for the second component in the map set I. 
Parameters are rescaled to normalize to unity in the Ql band 
(where the component dominates) in order to fix the degen- 
eracy between mixing parameter normalization and power 
spectrum amplitude. Although much lower than the CMB, 
this second component is clearly detected in the Q band 
(at about 10 standard deviations). Its mixing parameters 
are compatible with in V-bands, and are systematically 
negative and detected at about two standard deviations on 
average in each individual map of the W-bands except Wl. 
Mixing parameters for radiometers observing at the same 
frequency are compatible with a constant value, as expected 
for an astrophysical component. The amplitude of the com- 
ponent is at least 10 times larger in the Q band than in the 
V band (in Rayleigh Jeans (RJ) temperature units), and is 
larger in absolute value in the W band than in the V band. 
This behavior is not a priori excluded if the component orig- 



inate from residual galactic emission (mainly synchrotron 
and dust correlated emission) after galactic foreground re- 
moving operation by the WMAP team. In that case, mixing 
parameters of two different frequencies may have opposite 
signs depending on the sign of the error on the parameters 
of the template fit. 

Figure 1^ shows the map of the residual component ob- 
tained by Wiener filtering using estimated parameters. We 
can notice a bright structure close to the galactic center 
at the edges of the galactic mask which we identify as a 
residue from foreground subtraction in the Ophiuchus com- 
plex. Other residual structures associated to the North po- 
lar spur, the Gum nebula, the Orion-Endanus bubble and 
the Taurus region can also be identified, suggesting that the 
emission of these regions does not perfectly match the model 
used by WMAP team. 

We have found a significant correlation (« 40%) be- 
tween our second component and the WMAP synchrotron 
emission estimation at the Q band (using MEM, see Ben- 
nett et al. 2003b), which supports the Galactic origin of our 
second component. 

As can be seen from the Wiener map (and as ex- 
pected for galactic emission), the component looks very non- 
stationary over the sky. As the stationarity of components is 
assumed in our model, errors on mixing parameters proba- 
bly are somewhat underestimated but we do not expect this 
to change our interpretation. 

Figure |n] shows the estimated power spectra of the sec- 
ond component in both map sets I and II. For the estimation 
with the map set II, we have fixed the mixing parameters 
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Figure 5. Map of the residual component as "seen" in tlie Q-band, obtained by Wiener filtering using estimated parameters. 



of the residual component, as well as the amplitude of point 
source power spectrum to the results obtained in the map set 
I. This allows to compare the two power spectra estimates 
of a component with a fixed effective emission law. The es- 
timated power spectrum of the residual component in map 
set I is measured with high significance for multipoles be- 
tween ^ = 2 and I ~ 150 and form a plateau {l^c{t) « const) 
with an amplitude « 10—12 ^K^ in the Q-band. For higher 
multipoles, the beam transfer function of Q-band detectors 
strongly suppresses the signal, and the sensitivity is not suf- 
ficient to allow an accurate measurement of the component. 
Nevertheless, the last three bins show weak incompatibilities 
with at 1.5-2 sigma. Errors on those three parameters are 
strongly correlated (to « 50% due to the anti-correlation of 
those parameters with the amplitude of point sources emis- 
sion), so there is not clear evidence of detection at high 
spatial frequencies. 

The power of the residual component for angular scales 
larger than 1-2 degrees is clearly reduced at galactic latitude 
higher than 40 degrees. It remains marginally detectable 
even though it is reduced by a factor greater than 4 for 
I between about 10 and 30. 

We find that the power of the remaining galactic fore- 
ground contamination in WMAP maps (outside of Kp2 
mask) is about 1% of the CMB anisotropy variance at 
large angular scales {I < 100) at 41 GHz. This result is 
in agreement with the WMAP team estimate, based on cor- 
relation measurement using external foreground templates 
jBennett et alJl2003bl) . The power of the residual compo- 
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Figure 6. Estimated power spectrum of the second component 
from SMICA analysis in map set I (red ink) and comparison with 
the power spectrum estimated at galactic latitudes higher than 
40 degrees (map set II, blue ink). Power spectra amplitudes are 
rescaled on the Q-band. No positivity constraints are enforced for 
parameter estimates, making it easy to check their compatibility 
with 0. Cosmic variance, naturally included in error bar com- 
putation from the Fisher information matrix, has been removed 
afterwards from the error estimates. Estimated power spectra are 
re-binned for readability. Due to the lower sky coverage fraction 
for map set II (30%), the large angular scales are poorly con- 
strained. We then restrained the analysis to multipoles £ > 10 
only. 
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nent is less than 0.2 % of the CMB variance at 61 and 94 
GHz. 



4.4 Unresolved point source emissions 

In our analysis, the third component is rigidly constrained: 
as described above (sec. 14. it . the mixing parameters and the 
shape of the power spectrum are fixed to values expected for 
unresolved point sources emission, leaving only the overall 
amplitude to be determined by spectral matching. We find 
Cps(^) = (9.2 ± 5.0) X 10"^ ^iK^ at 41 GHz. This value 
is marginally compatible with the WMAP team estimate 
(seejiomatsu ct al. (2003) for a discussion for the expected 
contribution of unresolved point sources to the CMB power 
spectrum, f or different ma sks and fiux limits) and model 
predictions llArgiieso et aLl EoOSl: 15.5 ± 1.7 x 10"^-^'^- 

The large error bar on our estimate is due to the 
difficulty of discriminating components with similar spec- 
tral characteristics. In this experiment, the residual galactic 
component and the point source component have similar 
effective emission laws (mixing matrix elements) and both 
components dominate in the Q-band and can not be accu- 
rately estimated in V and W bands. So, some power at high 
spatial frequency can be exchanged between the second and 
the third components without modifying dramatically the 
likelihood of the model. 

We further investigate the distribution of the third com- 
ponent by performing two additional spectral matches on 
map set II. Neglecting residual galactic emissions (we have 
seen in previous section that residual galactic emission is 
clearly reduced at high galactic latitude) , we now fit a two- 
component model (as opposed to previous matches, build 
with three components) with two different constraint set. 
In a first match, all parameters related to point sources 
emission law and power spectra are fixed, except the am- 
plitude which is fitted. We find the amplitude of the power 
spectrum to be: Cps(^) = (11.3 ± 3.5) x 10"^ tiK^ at 41 
GHz, in better agreement with the value already reported 
of 15.5 ± 1.7 x 10~^ fiK^. In a second match, we relax all 
constraints on the shape of the power spectrum to obtain 
estimates of the power spectrum in every multipole bin. 
The resulting power spectrum of unresolved point sources 
is displayed in figure |7] Parameters are clearly incompat- 
ible with for most of the multipoles (assuming a priori 
a smooth power spectrum). A good compatibility with the 
model CpsW = 11.3 x 10"^-?^^ (41 GHz) is observed. If we 
re-adjust the model of the power spectrum to the estimated 
parameters using a simple we obtain a similar value for 
the amplitude. Nevertheless, a weak excess of power can be 
seen for the lowest multipoles {£ < 50) possibly due to a 
remaining galactic contamination at high galactic latitudes. 



4.5 Goodness of fit 

In this section, we briefiy examine the fit of the best models 
to the data across frequency bins by plotting the weighted 
mismatch nqD{Rq, Rq{9)) against the bin index q. If the 
model of independent stationary component holds, the ex- 
pected value of this measure of mismatch is given, for riq 
large enough, by eq. (1151 . In our plots, however, we use 
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Figure 7. Power spectrum of unresolved point sources emission 
estimated with SMICA assuming two components in the map 
set II (red). Mixing parameters of point source emission have 
been fixed for the analysis. No constrains of positivity is given 
for parameter estimates allowing to check easily their compati- 
bility with 0. Cosmic variance has been removed afterwards from 
the error estimates. The black curve corresponds to the model of 
point source power spectrum Cps (^) = const with an amplitude 
estimated from the semi-blind approach of 11.3 X 10~^iJ,K^ (see 
text for details). 



Monte-Carlo simulations to have an estimate of the distri- 
bution of the mismatch valid even in the non asymptotic 
regime. We report goodness of fit of adjusted models in 
two cases: a one-component model and a three-component 
model. 

FigurelHlshows that the one-component model is clearly 
incompatible with the data for £ < 60. It also shows that the 
fit obtained with three components is very satisfactory for 
most of the multipoles although small discrepancies remain, 
in particular for £ between 10 and 25. We interpret this as 
a clear indication that residual galactic emission are present 
in the data at low spatial frequency. 

We can investigate the origin of the remaining discrep- 
ancies by looking at the spectral mismatch between pairs of 
detectors. Figures |^ shows the goodness of fit between the 
individual pair of detectors (Q1,Q2), assuming either one or 
three components. The misfit for the one-component model 
is now even more obvious (as compared to the global mis- 
match involving all detectors) , as expected since the Q band 
is where the galactic residual dominates. 

As the second component is detected in the W channel 
as well, we also looked at the fit of one-component and three- 
component models for all W detectors together and for the 
set of W2, W3 and W4 (discarding Wl). For both figures, 
the goodness of fit is significantly better for three compo- 
nents than for a unique component for low multipoles. 

A significant excess remains, which is probably due to 
fine departure of the residual galactic emission from the sin- 
gle template assumption. Figure 1101 however, rather point 
to remaining systematics in the Wl channels since the mis- 
match is significantly reduced when excluding the Wl de- 
tector. 

This incompatibility, however, is small enough so that 
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Figure 8. Spectral mismatch for map set I, with three compo- 
nents (red) and 1 component (blue). The green curves are the 
boundaries of the 68% goodness of fit interval estimated using 
simulations. For three components, the fit is very satisfactory for 
most of the spatial frequencies (see text). 
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Figure 9. Spectral mismatch for the detector pair (Q1,Q2) only, 
with three components (red) and one component (blue). The 
green curves are the boundaries of the 68% interval estimated 
using simulations. Note the logarithmic scale and thus the very 
significant reduction of the mismatch when three components are 
considered. 



the CMB power spectrum is not affected because the errors 
for the lowest multipoles are dominated by cosmic variance. 



5 CONCLUSIONS 

We have performed a blind multi-component analysis of 
the first year WMAP data. We used the eight foreground- 
cleaned high frequency maps (Ql to W4) provided by the 
WMAP collaboration. Our analysis uses SMICA, a maxi- 
mum likelihood spectral matching method. Three significant 
astrophysical components 

• the CMB anisotropics 

• a residual galactic emission 

• unresolved point source emissions 
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Figure 10. Spectral mismatch with three components (red) and 
one component (blue). Top: mismatch for all four W band detec- 
tors; bottom; for three last W band detectors (Wl excluded). 



are jointly characterized in the foreground cleaned maps 
after masking the strongest point sources and the galac- 
tic plane using Kp2 mask. No significant thermal SZ 
emission is found with the present study, in a gree- 
ment wit h the WMAP team (iBennett et al.l l2003b|) and 
Huffcnbcr ger. Seliak. Makarovl ll2004f) . 

Blind analysis allows us to estimate the power spectrum 
of CMB anisotropics as well as their amplitude across fre- 
quency channels confirming their cosmological nature. Our 
power spectrum estimate is in excellent agreement with the 
WMAP team estimate. We show that the measured emission 
law of CMB anisotropies at WMAP frequencies is compati- 
ble with the derivative of a blackbody, as expected for tem- 
perature fiuctuations. The statistical errors on parameters 
related to the amplitude of the anisotropies are about 0.3%, 
and are smaller than calibration errors (0.5%) on dipole 
modulation. Conversely, if CMB anisotropies are assumed 
to be pure temperature fluctuations, then the estimation of 
CMB amplitude across channels provides a relative calibra- 
tion across frequency bands at a better precision than dipole 
calibration. 

The second estimated component, corresponding to a 
weak residual galactic emission, is mainly concentrated in 
Q-band maps outside of the Kp2 mask. We believe that 
this component results mainly from spatial variations of 
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the difference between Haslam map (used as a synchrotron 
template for the subtraction) and synchrotron emission at 
WMAP frequencies. This component is weak compared to 
CMB anisotropics. The estimated power spectrum is about 
£{£ + 1)C{£)/2tt f» 10 - 12fiK^ for £ < 100 in Q band and 
less than 2^,K^ in V and W band. Those estimates are com- 
patible with WMAP team estimates of foreground contam- 
ination. Finally, much of the power from this residual large 
scale component disappears for galactic latitude higher than 
40 degrees. 

The third component corresponds to residual point 
sources emission. By fixing the parameters related to the 
emission law of point sources (we assume {u/vq)~^) and to 
the power spectrum (flat power spectrum), we estimate the 
amplitude of unresolved point source power spectrum. We 
find at high galactic latitude: C{C) = (11.3±3.7) x IQ'^ ^.K^ 
at 41 GHz compatible with the WMAP team estimation. 
We also provide an estimate of the power spectrum of point 
sources at high galactic latitude. 

The goodness of fit of the 3-component model is ex- 
cellent, except for a very small discrepancy around ^ ~ 20 
which is likely to be attributed to a systematic in Wl. This 
inconsistency has no impact on the CMB power spectrum 
estimate (because errors are dominated by cosmic variance), 
but could affect the estimates of weak components in data. 
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